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We theoretically study the dissipative dynamics of a quantum emitter placed near the planar 
surface of a metal supporting surface plasmon excitations. The emitter-metal coupling regime can be 
tuned by varying some control parameters such as the qubit-surface separation and/or the detuning 
between characteristic frequencies. By using a Green's function approach jointly with a time- 
convolutionless master equation, we analyze the non-Markovian dissipative features on the qubit 
time evolution in two cases of interest: i) an undriven qubit initially prepared in its excited state 
and a) the evolution towards a steady-state for a system driven by a laser field. For weak to moderate 
qubit-metal coupling strength, and on timescales large compared to the surface plasmon oscillation 
time, a Markovian approximation for the master equation results to be adequate to describe the 
qubit main optical properties: surface enhancements of rate emission, optical spectra and time- 
dependent photon-photon correlation functions. The qubit decay shows a crossover passing from 
being purely dissipative for small qubit-surface distances to plasmon emission for larger separations. 



I. INTRODUCTION 

Surface plasmons (SP) on metals, a topic exten- 
sively studied from many years ago^, has recently re- 
ceived a strongly renewed attention due to significant 
advances in new experimental capabilities and numerical 
development^^. Great attention has been focused on 
the emerging field of quantum plasmonic with the goal 
of making devices for quantum information processing^ 
as single-photon transistoi!^ or laser^. As a requisite for 
this goal, a lot of effort has been devoted to get coherent 
coupling between plasmons and a quantum emitter made 
of a solid state qubit (SSQ) as, for instance, a quantum 
dot, a single NV center or a single molecule among oth- 
ers. 

Strong coupling signatures of SSQ and SP have been 
experimentally observed both in spectroscopic as well 
as in time-resolved studies. In spectroscopy, the anti- 
crossing between exciton and plasmon features in opti- 
cal spectra have already been reported in systems such 
as organic semiconductors^, organic molecules placed in 
subwavelength hole arrays^, metallic nanowires^, hybrid 
metal-semiconductor nanostructures^ and even in car- 
bon nanotubes^^. These anticrossings have been claimed 
to be a manifestation of strong coupling between SSQ and 
SP. On the other hand, ultrafast time-resolved signatures 
of strong coupling in SSQ-SP systems have also been re- 
ported. An enhancement of several orders of magnitude 
for the spontaneous emission rate in a time-resolved pho- 
toluminescence measurement on a InGaN heterostruc- 
ture close to a silver thin layer has been reporte d^^ ' ^^' . 
Additionally, recent experiments which operate simulta- 
neously with both Raman and fluorescence signals com- 
ing from a single molecule in very close proximity to a 
metal surface have allowed the indirect measurement of 
ultrafast 25 fs) dynamical features in such SSQ-SP 
systenP^. 



On the theoretical side, some progress has been made 
to understand SSQ-SP coupling in different geometries 
using different approaches. The first attempts were de- 
voted to computing the spontaneous emission rate en- 
hancement of an atom near an absorbing surface as given 
by the atom self-energy in a near field limiPl. More re- 
cently, a hydrodynamic model^^ has been used to study 
a single molecule coupled to metallic nanoparticles. A 
transfer matrix method has also been used for simulat- 
ing attenuated-reflection experiments^^. However, the 
experimental setup which has risen the highest inter- 
est has bee n the quantum emitter coupled to a metal- 
lic nanowir^^^^, where the generation of a single optical 
plasmon can be achieved. Several theoretical studies on 
this system have considered the full quantum behavior of 
plasmon modes^^ In particular, some attention has 
been devoted^ ^ to non-Markovian effects that can be im- 
portant in the SSQ-nanowire system because the spectral 
density J{uj) (carefully discussed in the present work) is 
highly structured due to a divergence at the edge of the 
SP density of states. 

An open quantum system strongly coupled to a 
reservoir displays a complex dynamics which, in gen- 
eral, requires a description beyond simple Markovian 
theories^^ In order to clarify the relevance of non- 
Markovian effects in SSQ-SP systems, we concentrate in a 
quantum emitter close to a planar surface of a dissipative 
metal, a system conceptually simpler than wires because 
it only has a single band of plasmons with a density of 
states having a singularity at a frequency ujsp- We study 
the properties of the light emitted by the system depicted 
in figure [l] a SSQ close to the planar metallic surface 
which supports a plasmon field as well as some dissipa- 
tion mechanism. Strong SSQ-SP coupling could be ex- 
pected when the qubit-surface distance is small compared 
with a typical length scale as, for instance, the wave- 
length of the emitted light. In order to understand the 
fundamental mechanisms of SSQ-SP strong coupling, we 
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restrict ourselves to consider just a single quantum emit- 
ter. However, collective effects of many emitters coupled 
to the same plasmon field have been recently proposed^ 
as responsible for the detection of the Rabi vacuum split- 
ting in these systems. 




FIG. 1: Schematic view of the system: a SSQ, with charac- 
teristic frequency cjq, is placed at a distance z of an infinite 
planar metallic surface. 

We start by analyzing non-Markovian features in the 
SSQ-SP optical features by using a time-convolutionless 
approaclF^. We show that SSQ-SP dissipative evolu- 
tion is determined by the rapidly varying structure of 
the reservoir spectral function at a frequency close to 
cjgp of the order of a few eV in a normal metal. There- 
fore, the timescale for non-Markovian effects reduces, at 
most, to a few hundreds times eV~^^ i.e. typical times 
under a picosecond. Consequently, we conclude that the 
observation of non-Markovian signatures in SSQ-SP sys- 
tems made with normal metals will indeed demand ex- 
perimental set ups at the edge front of present state-of- 
the-art ultrafast technology. Beyond this short timescale, 
rate emissions stay constant and a Markovian approach 
becomes adequate for computing population dynamics, 
optical spectrum and second order coherence function. 
Here, both a Markovian and a non-Markovian analysis 
are applied to a SSQ-SP system under two different exci- 
tation schemes: firstly, we analyze the spontaneous emis- 
sion of a SSQ initially prepared in the excited state. Sec- 
ondly, the time evolution of a SSQ initially in its ground 
state and driven by means of a coherent laser field up to 
a stationary state is studied. 

We take h = 1 along this paper which is organized as 
follows: in Section [ill we introduce the Green's tensor of 
the layered system and study its main properties, in par- 



ticular the spectral density function. In Section III the 
time-convolutionless method is briefly reviewed and the 
non-Markovian effects on the SSQ-SP system dynamics 
are considered. In Section [TVl we use the Markovian limit 
to calculate the optical properties of the system. Finally, 
in Section [Vl we summarize our results and draw some 
conclusions. 



We study the system depicted in figure [Tj a SSQ in the 
upper-half space is embedded within a dielectric matrix 
with a dielectric function that can be taken as real and 
constant, ei, in the range of frequencies of interest. In 
the lower half-space, z < 0, a dissipative metal is char- 
acterized by a complex dielectric function 62 (cj) that we 
take in a renormalized Drude approximation: 



e2(^) 



(2) 



600 is the high-frequency limit of the metal dielectric func- 
tion. Up is the bulk plasmon frequency and jp is the Lan- 
dau damping constant. 




FIG. 2: Spectral density J{uj)/jo for a detuning A = 
Cosp — Coo = 0.1 plotted for different qub it-surface separa- 
tions: z = Upz/c = 0.01 (blue), 0.32(purple), 0.64(yellow) and 
1.42 (green). Notice that as far as the SSQ approaches the 
surface, the spectral density increases in roughly 6 orders of 
magnitude. 



The Green's tensor for this layered geometry has two 
contributions: 

G(r, r', uj) = Go(r, r', oj) + Gnir, r', oj) 

where the first term is the free-space solution given 
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II. GREEN'S TENSOR AND SPECTRAL 
DENSITY 

Electrodynamics of a dissipative medium is described 
by the Green's tensor G(r, r',co') which satisfies the 
Maxwell equation: 
. ,2 



(V X Vx) 



re(r,cj) 



G(r,r^c^) = I(5(r-rO. (1) 



with q = uj/c^ (j){x) = x-\-ix'^ — x^, x(^) = x-\-3ix'^ — 3x^, 
I is the 3x3 identity matrix and R (g) R represents the 
dyadic product between the vectors R = r — r^ The sec- 
ond term in Eq.([3|, Gi^(r, r', cj), is the reflection Green's 
tensor with on-site nonzero components^ 
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where tii{q^uj) = q^ — ei{(jj){(jj/cY for z = 1,2. 

All the parameters of the absorbing medium relevant 
to the SSQ dissipative dynamics appear in the Green's 
tensor G(rQ, rg, cj), where yq denotes the SSQ location. 
The action of the absorbing medium on a SSQ with dipole 
moment p, is completely described by the spectral den- 
sity 



J{uj) 
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Im[G(YQ,YQ,uj)] 



(6) 



which is related with the qubit-environment coupling 
g{uo) and the density of states of the environment p{(jj) 
by means of J{uj) = g'^ {uj) p{uj) . In order to compute 
the spectral function of a representative system, we use 
parameters for silver in the range of frequencies of inter- 
est where ujp = 3.76eV^ e^o = 9.6 and 7^ = O.OScjJ^, 
while for the dielectric constant at the upper-half-space 
we take ei = 5. S ince the densit y of states has a singular- 
ity at ujsp = ^p\/^oo/{^i + ^oo) = O.Slujp non-Markovian 
effects associated with the structured reservoir can be 
expected to occur around that inverse frequency. Thus, 
we consider a SSQ with a dipole oriented along the z- 
direction and an energy splitting uq slightly detuned with 
respect to the singularity at the SP edge uJsp. In particu- 
lar, we calculate J{uj) for A = Osp — ujq = 0.1 where the 
frequencies have been normalized to the bulk plasmon 
frequency COsp = uJsp/uJp and Coq = uJo/ujp. 

In figure [2] we plot, for different values of the qubit- 
surface renormalized separation z = zup/c^ the spectral 
density in units of the spontaneous decay rate of the SSQ 
in free space 70 = uJqp'^ /Sttcoc^ . Two main results can 
be identified: 

i) A strong reduction of J{uj) when the qubit gets far- 
ther from the surface. This is a consequence of the expo- 
nential reduction of the coupling, as a function of z^ as 
indicated in Eq.Q. 

ii) For small separations J{uj) is highly structured pre- 
senting a strong peak close to the frequency ujgp (vertical 
line in Fig.Q). This is a consequence of the singularity 
of the density of SP states at small detunings A. For 
increasing separation, J{uj) becomes much smoother and 
a reduced maximum separates from uj sp- 



ill, TIME CONVOLUTIONLESS METHOD AND 
NON-MARKOVIAN EFFECTS 



The spectral density for the SSQ-SP system computed 
within a purely classical scheme, can be used within a 
quantum framework describing the dissipative dynamics 
of an open quantum system. As it is well known, when 
the time correlation between the system and the environ- 
ment decay much faster than the characteristic inverse 
dissipation rate, memory effects can be neglected in the 
so called Markovian approximation, and the observables 
of the system are given by analytical expressions. How- 
ever, this timescale does not represent the only one rele- 
vant to determine the system's evolution. When the envi- 
ronment correlation time is longer than the inverse rate of 
the system-bath coupling, new physics can arise at very 
short times. A SSQ in close proximity to a dissipative 
metal surface supporting SP modes sees a highly struc- 
tured reservoir to which it might be strongly coupled. 
Thus, one can expect non-Markovian effects to be signif- 
icant in the qubit time evolution. Many works^^ 29lEo]^ 
have been devoted to treat this problem at different lev- 
els of precision and sophistication. Here, we chose to 
work within a time-convolutionless (TCL) frameworlP^ 
to capture non-Markovian effects to the lowest order in 
the SSQ-SP coupling strength. This method has already 
been applied to consider the spontaneous decay of a two- 
level system coupled to a general structured reservoir^. 
For SSQ-SP systems, the strong variation of J{uj) occur- 
ring for frequencies close to ujsp^ implies that dynamical 
features in time scales from femtoseconds to picoseconds 
are expected. 



A. Time-convolutionless method 



What is of interest for us of the TCL method can be 
sketched as follows^^. It consists in transforming the typ- 
ical non-Markovian integro- differential equation for the 
reduced density matrix into a local in time evolution 
equation by making use of a power expansion technique of 
the Nakajima-Zwanzig type. As a result, a master equa- 
tion for a qubit is obtained with time dependent decay 
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rate 7(t) and Lamb shift S{t): 
dp{t) i 



dt 



+ i^(2a-p(t)a+ - cj^cj-pd) - p{t)cj^cj-) (7) 

Time-dependent rates can be calculated within a pertur- 
bative expansion. In order to calculate them, a first step 
is to Fourier transform the spectral density: 



fit) = J dujJ{uj)e 



(8) 



The lowest order non-Markovian effects, i.e. the so-called 
post-Markovian behavior, are contained in the second or- 
der contributions to ■~f{t) and S{t) given by: 



1 /■* 

72 W = dti^f{t-h) 
1 /■* 

S2{t) = ^J^ dh'^f{t-h). 



(9) 
(10) 



where 3? and ^ denote real and imaginary parts, respec- 
tively. 



B. SSQ spontaneous decay 



resonance with the top SP energy, e.g. A = 0.1, the os- 
cillations slow down, the envelope of the oscillatory decay 
rate becomes smaller and the negative parts of the decay 
rate 72 (^) tend to vanish producing only few oscillations 
before the spontaneous decay becomes almost exponen- 
tial. For further smaller detunings, e.g. A = 0.01, the 
SSQ sees an even more structured reservoir with a de- 
cay rate modifying completely its behavior: its value in- 
creases considerably and it just oscillates slightly around 
a large positive value, producing a monotonous decay of 
the SSQ excited state population. This last result indi- 
cates that the second-order TCL method is approaching 
its limit of validity. Physically, this behavior is a conse- 
quence of the fast transfer of the SSQ energy to the SP 
field, an energy which is irreversibly lost and the quan- 
tum emitter ends up in its ground state. 
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We start by considering the situation where an un- 
driven SSQ is prepared in its excited state from which 
it decays emitting a photon to the vacuum or to the SP 
field. There are two possible situations depending on the 
sign of the detuning A between SP COsp and the SSQ ooq 
renormalized frequencies. The time evolution of the ex- 
cited state population, is given hy '^^ 

with ni(0) = 1 and the decay rate obtained from Eq.([9|. 

At this stage, we want to analyze the importance of 
memory effects. Therefore, in the calculations reported 
in this subsection we do not include the free space part of 
the Green's tensor, Eq.([3|, which involves a much slower 
dynamics than the one associated to the reflection con- 
tribution Gr^ as given in Eq.Q. 

Figure [3] shows ni{t) for different positive detunings, 
i.e. when the SSQ is resonant with the continuum stripe 
of SP modes(0 < uj < oJsp)- In order to have a highly 
structured reservoir, we have taken a small qubit-surface 
separation, z = 0.055 (very close to the blue line spectral 
density in figure |2|. For large detuning A = 0.5, 72 (^) 
oscillates around a constant (Markovian) value. At some 
time intervals, 72 (^) takes on negative values, a fact that 
tends to slow down the decay of the excited state pop- 
ulation. Physically, this behavior can be understood as 
due to the back-action of the reservoir on the SSQ re- 
exciting it. When the SSQ splitting energy gets closer in 
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FIG. 3: Decay rate (a) and population of the excited state 
(b) of a SSQ located at z = 0.055 from the planar surface. 
Different lines correspond to different detunings from the SP 
frequency: A = cj^p— cjq = 0.5(solid blue), 0.1(dashed purple) 
and 0.01 (dotted yellow). 

On the other hand, a physically different situation oc- 
curs for negative detunings, i.e. when the SSQ energy 
falls within the gap where no SP states exist. Figure 
|4] is similar to figure |3] with the same z and detunings 
with just a change of sign with respect to the ones in 
figure [3| When uq is far above the edge of the SP dis- 
persion, the SSQ basically remains in its excited state 
as the spectral density for these energies is practically 
zero, so there are not accessible SP modes to which de- 
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comes constant just at the value it takes in a Markovian 
description as discussed in the fohowing Section IV. 



C. Coherently driven SSQ 

After having studied the effect of the structured reser- 
voir on the SSQ spontaneous emission, now we turn 
our attention to the case where the system is coherently 
driven by a laser field. The SSQ emits and absorbs pho- 
tons simultaneously. The system can achieve a stationary 
state in which light absorbed from the laser ends being 
transferred to plasmons. The laser can be treated as 
a classical field included in the, local in time, coherent 
part of the master equation ^ through the Hamiltonian 
i^ias = l^(cr+e^^i-^ + e-^^i-V-)/2. The Rabi frequency 
Q measures the strength of the coupling of the SSQ dipole 
with the laser field. The time evolution of the SSQ ex- 
cited state population is obtained from the solution of 
Eq.Q represented in a rotating frame at the laser fre- 
quency cjias- ^i(^) is plotted in figure [sj for the case of 
perfect resonance between the laser and the SSQ. We as- 
sume that 72 (^) remains unaffected by the laser field, so 
that the main effect of the laser is to bring the system 
to a stationary state in a timescale which is similar to 
that of the spontaneous decay discussed in the previous 
subsection. 



FIG. 4: Decay rate and population of the excited state for 
a SSQ placed at z = 0.055 from the planar surface. The 
different lines correspond to different detunings: A = Cosp — 
Coo = — 0.5(solid blue), -0.1(dashed purple) and -0.01(dotted 
yellow). Inset in (a) corresponds to a zoom of the decay rate 
for the first two values of A, while inset in (b) corresponds 
to a zoom to population of the excited state at short times 
for A = —0.5. Including the effect of free space emission (70) 
produces a decay of ni in a time scale much larger than that 
of the figure. 



cay into. Nevertheless, as shown in the inset of figure 
[3j one may observe some non-Markovian oscillations for 
very short times. When the SSQ energy is tuned closer 
in resonance with COgp^ e.g. A = 0.1, a very interesting 
phenomenon occurs: the emitter undergoes the so-called 
fractional decay in which the population tends to a finite, 
non-zero, value at long times. Including the effect of free 
space emission (Gq) produces a decay of rii in a time 
scale (70""^) much larger than that of the figure. As the 
emission frequency is further scanned closer to the band 
edge, the behavior changes again dramatically: the decay 
rate, instead of oscillating around zero, oscillates slightly 
around a positive value, which results into an irreversible 
exponential decay. 

A very important result must be drawn from all these 
results: the timescale of these non-Markovian effects is a 
few hundred times oj~^. For normal metals, this means 
times below Ips. Beyond that short timescale, 72 be- 



IV. OPTICAL PROPERTIES IN THE MARKOV 
APPROXIMATION 

The results of the previous section show that a SSQ 
presents significant non-Markovian effects in a timescale 
a couple of orders of magnitude larger than UJ~^ . Here- 
after, we concentrate in the usual case of having a reso- 
lution in time larger than a picosecond. Then, the sys- 
tem can be described by a Markovian dynamics given 
by a Master equation like Eq. ^ but now with a 
Lamb-shift S = S2{t ^ 00) and a constant decay rate 
7 = 72(^ ^ 00) including both the free space and the re- 
ffection contributions to the dissipative dynamics. Since 
the only effect of the Lamb-shift is a constant energy 
shift, from now on we do not pay attention to it. 

A. Decay rate 

The decay rate in the Markovian regime coincides with 
the long time limit of 72 (t), 7 = 72(^ ^ 00), allowing to 
identify 7 as simply the spectral function at the SSQ 
frequency: 

7 = 2^J(c^o) = • /m[G(rg,rg,cjo)] • P (12) 

where the two terms corresponding to the free space (Go) 
and the reflection part (G^^) of the dissipative dynamics 
are included in G. 
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FIG. 5: Decay rate and excited population of a SSQ placed 
at z = 0.2 with energy detuning A = Cosp — Cuo = 0.1 and 
coherently driven by a laser in resonance with the SSQ. The 
different lines correspond to different laser intensities Q/jo — 
lOO(blue), 200(purple) and 500(yellow) 



The SSQ decay rate to the SP reservoir of the metallic 
surface, 7, is shown in figure |6j in a parameter space 
{z^ujo}^ where lighter blues correspond to high decay 
rates with a variation of four orders of magnitude be- 
tween the highest and the lowest values. In order to 
discuss these results, it is better to plot 7 vs. the SSQ- 
interface distance z (in logarithmic scales) for different 
SSQ energies as depicted in figure [7| It is worth notic- 
ing two important features: First, at a large ujq value 
the assisted decay rate is smaller than the vacuum one 
for a certain range of distances, due to the fact that the 
reflected part of the Green's tensor is interfering destruc- 
tively with the direct one. This effect is evident when the 
SSQ frequency approaches Usp while it moves to larger 
separations and it weakens, when cuq is far from the 
SP band edge. 

Second feature is even more important. When the 
SSQ-surface distance varies, the decay rate suffers a tran- 
sition from a 1/2;^ behavior to a much slower decrease. 
In order to understand the physics behind this behav- 
ior, we analyze the much simpler case ei = 600 = 1 and 
7p = cjp/500, i.e. a rate for losses one order of mag- 
nitude smaller than the one we have used in previous 
cases. Large dots in figure [S] depict the decay rate as 




0.60 0.65 0.70 0.75 
UJQ 



0.80 



0.85 



FIG. 6: Surface assisted decay rate log^Q (7/70) given by Eq. 
( 12 ) for a region of the z, Cuo parameter space. The red dotted 
gridline marks the cjsp frequency. Lighter blues correspond to 
high values. Between the highest and the lowest values there 
are four orders of magnitude. 
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FIG. 7: Logarithmic relative decay rate logiQ(7/7o) calcu- 
lated with Eq.( |12| ) as a function of \ogiQ{z). The three curves 
correspond to different SSQ energies Coq = 0.24(solid blue), 
0.4(dashed purple) and 0.76 (yellow dotted). 



calculated with the full Green's tensor for two different 
SSQ frequencies. At very short distances the decay rate 
is highly enhanced with respect to 70 and it shows a 1 /z^ 
dependence, which can be obtained (dotted lines in the 
figure) by means of a model which only contains non- 
radiative processes as the creation of electron-hole pairs 
in the metallic medium. A crossover to a different be- 
havior occurs at a critical distance, which depends on 
the SSQ-SP detuning A. Beyond this critical distance, a 
single plasmon pole approximation (dashed lines in the 
figure) in the Green tensor in Eq.( 12 ) is able to reproduce 
the numerical result with the complete Green's tensor. 
In other words, for separations beyond the crossover, the 
SSQ decay just produces the emission of surface plasmons 
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while other metahic losses become negligible. In the case 
of the SSQ embedded in a dielectric or a metal with very 
large losses this crossover can be hindered by other phys- 
ical effects such as those coming from local dissipative 
circulating currents^^. As the SP channel contribution 
increases when the SSQ energy gets closer to the plas- 
mon band edge, this crossover effect can be exploited in 
designing coherent plasmonic devices^ ^ . 



B. Spectrum and correlation functions 

Hereafter, we consider the case in which the system is 
coherently driven by a laser so that the system reaches a 
stationary state with partial occupation of the two levels 
of the SSQ. The master equation in the rotating frame 
at the laser frequency, cJias, is: 

^ = z{Ac.[p(t),a+a-] + + a']} 

+ ^(2a-p(t)a+ - a^a-p{t) - p{t)(j^(J-) (13) 

with Ao; — uoq — cjias- The problem reduces to the study 
of the SSQ resonance fluorescence^ near the planar sur- 
face of a dissipative metal. The vacuum resonance flu- 
orescence has been widely studied in the literature in 
the case of resonant excitation for which an analytical 
solution for the population, spectrum and second order 
coherence function exists. Here, we extend such anal- 
ysis to the non-resonant case and pay special attention 
to the effect of the SP reservoir. We present our analy- 
sis in the way the experiments can be performed either 
by tuning the laser resonantly with the SSQ energy and 
varying the laser intensity or by fixing the laser intensity 
and scanning the laser frequency. Figure [9] explores the 
former alternative whereas in figure [lO] we consider the 
later one. 



From the master equation (13) one may derive the 
equations of motion for the expectation values (cr+(t)), 
{a~{t)) and {a^{t)) arriving to the well-known optical 
Bloch equations(OBE). The steady-state solution for the 
excited state population is: 



^2 



(72+4Aa;2 + 2cj2)- 



(14) 



In general, the OBE must be solved numerically in order 
to get the population dynamics (ni(t)), except for the 
resonant case (Acj = 0) for which an analytical solution 
exists: 



(niit)) 



7 + 2f22 

I cos{m) 



37 
4R 



sm{m) 



(15) 



where R = y^Q'^ — 7^/16, labeled as Rabi splitting at 
resonanc^', is the parameter characterizing the strength 



of the effective coupling. There is a threshold for the laser 
intensity at 1] = 7/4. For Q below this threshold, the so- 
lutions are monotonically decaying functions of time so 
that the system is said to be in the weak coupling (WC) 
regime. Above that threshold, the populations exhibit 
oscillations, and the system is said to be in the strong 
coupling (SC) regime. In figure |9|a) we plot the real 
part of R in the parameter space {z^CIq}. The bluest 
region corresponds to 3?(i?) = 0, which means that the 
Rabi splitting at resonance is purely imaginary and con- 
sequently the system is in the WC regime. For the regions 
in which blue becomes lighter, the values correspond to 
positive and higher values of J?(i?). 

In order to clarify these results, we show the popu- 
lation dynamics in figure 9jb) for three different points 
highlighted in part (a) of the same figure: the green curve 
corresponds to a configuration where the laser is weakly 
coupled to the system, so no oscillations are observed in 
the population. The red point corresponds to the region 
of transition from WC to SC where just one clear os- 
cillation occurs before practically arriving to the steady 
state. Finally the blue point corresponds to a configura- 
tion where the laser is strongly coupled to the SSQ and 
several oscillations are observed before the steady state 
is achieved. 

Another experimental alternative is to keep Q constant 
and vary the laser frequency as it is plotted in figure 10 
In this case, the laser is out of resonance and the Rabi 
splitting must be redefined as^^: 



Ra = 



(16) 



Figure 10 ^a) shows 3?i?A in the parameter space {^, AcD}. 
The population dynamics is shown in panel (b) of the 
same figure, where one can observe the transition from 
strong coupling (oscillations, solid blue) to weak coupling 
(monotonous in time, green dotted) for non-resonant ex- 
citation of the SSQ. 

A clear manifestation of the transition from WC to SC 
appears in the optical spectrum at the stationary regime. 
It can be calculated through the Wiener-Khintchine for- 
mula: 



1, 



S{uo)) = -3? / e*^^(a+(t)a-(t + r))dr. (17) 
7^" Jo 

The calculation of the two-time correlator in Eq. ( [iTt 
requires the use of the Quantum Regression Theoren^ 
by using the steady state populations as initial values 
for the second time dynamics. In the resonance fluores- 
cence problem there are always two main contributions 
to the spectra: the Rayleigh scattering coherent part and 
the one coming from the incoherent scattering. The for- 
mer contribution is just a delta function at ujq that we 
ignore in our results. We are mainly interested in the 
contribution coming from the inelastic scattering which 
is shown in figures |9|c) and [lOj^c) for the resonant and 
non-resonant cases, respectively. As it occurred with the 
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population, under resonant excitation the spectrum ad- 
mits an analytical expression: 

7/4 



S{co) oc ( 



37/I6 



37/16 



■ Wo 



(18) 



In the WC regime (green point) Q < 7/4, the light emit- 
ted simply produces a Lorentzian curve peaked about uq 
with linewidth 7/2. For the intermediate regime (red 
point), on top of the Lorentzian peaked at the qubit 
frequency, some satellites start to appear at the laser 
Rabi frequency ±1^. For a strong-driving field situa- 
tion Q > 7/4 these two sidebands appear at frequen- 
cies LjJ = cuq ^ Q. For the non-resonant case, the thresh- 
old changes but the behavior remains qualitatively unaf- 
fected: even though the dressed state structure is slightly 
modified by the detuning, at the end, a triplet is ob- 
tained in the resonant case. The existence of this Mol- 
low's triplet is a manifestation of the SC of the laser to 
the SSQ-SP system. 

Another magnitude of experimental interest is the sec- 
ond order coherence function: 



G(2)(t,t + r) 



G(i)(0)G(i)(t) 
with correlation functions 



(19) 



G(2) {t, t + T)= ((t(+) {t)a^+^ {t + T)a^-\t + T)a^-\t)) 
GW{t,t + T)= ((T(+)(t)cT(-)(t + T)). (20) 

We evaluate these magnitudes at the stationary state. In 
the resonant case, the second order coherence function 
can be analytically expressed as: 



-37r/4 



cos{Rr) 



37 
4R 



sm{Rr) . (21) 



It clearly exhibits photon anti-bunching: ^^(0) = 0. Fig- 
ure ^d) shows ^f^^^ for zero detuning for the three dif- 
ferent points considered above for the other magnitudes. 
Apart from the antibunching, the case of SC shows a re- 
markable oscillatory behavior. Once more, qualitatively 
similar results are obtained with laser-SSQ detuning as 
shown in figure [lOj^d) . 

The main consequence to be drawn from figures [9] and 
10 is that by pumping the SSQ with a tunable laser, 
and measuring spectra and second order correlation func- 
tions, one can extract information about the SSQ cou- 
pling to the surface plasmon of the dissipative metal. 



V. SUMMARY 

In this work we have studied the properties of the cou- 
pling of light with a SSQ, embedded in a dielectric, in the 
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FIG. 8: Logarithmic relative decay rate log;Lo(7/7o) calcu- 
lated with Eq.(|12[) as a function of \ogiQ{z) with ei = eoo = 1 
and 7p = cjp/SOO. Two different SSQ energies are consid- 
ered: cjo = 0.2 (blue circles, far from the SP edge) and 0.6 
(red circles, close to the SP edge). Dashed lines denote the 
decay rates as predicted by a single pole approximation (SP 
contribution). Dotted lines correspond to metal losses as the 
emission of electron-hole pairs. 



presence of a SP field supported in the interface between 
this dielectric matrix and a dissipative metal. Using a 
time-convolutionless approach, we provide a theoretical 
description of the non-Markovian features for this kind of 
systems and discuss its relevance in possible observations. 
In a spontaneous decay situation, different behaviors oc- 
cur depending on both the sign and the absolute value 
of the SSQ-SP detuning: from a monotonous (almost ex- 
ponential) decay for very small detunings, to population 
oscillations due to reabsorptions in the case of positive 
detuning. Even fractional decays can be observed, when 
negative detunings are present and the SSQ energy is not 
too close to the SP edge band. 

In experimental situations, non-Markovian features 
can be hard to detect due to practical difficulties in get- 
ting the adequate time-resolution. Therefore, we have 
also considered a Markov approximation to study the 
electrodynamics of the SSQ coupled to a reservoir of SP 
modes. The whole information of the planar metallic 
surface is embedded in the decay rate constant, which 
depends on both the SSQ frequency and distance to the 
surface. The excitation of the system by a laser allows 
the existence of a steady state as well as the analysis 
of different measurable properties of the SSQ-SP system 
as, for instance, surface enhancements of rate emission, 
optical spectra and time-dependent photon-photon corre- 
lation functions. Our main result is that the qubit decay 
shows a crossover passing from being purely dissipative 
for small qubit-surface distances to plasmon emission for 
larger separations. As the SP emission channel increases 
when the SSQ energy gets closer to the plasmon band 
edge, this crossover effect can be exploited in designing 
coherent plasmonic devices. Our next task, beyond the 
scope of the present work, is to treat the plasmonic part 
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FIG. 9: Optical properties of the SSQ-SP system. Panel 
(a) corresponds to the value of ^(R) in the parameter space 
{z, Clo} in order to distinguish the strong and weak coupling 
regions. Panel (b) and (d) show the dynamics of the ex- 
cited state population and the two-photon correlation func- 
tion {g'''^\t)) respectively for the three points plotted in panel 
(a). Panel (c) shows the qubit luminescence spectra for those 
three particular cases. 



FIG. 10: Optical properties of the SSQ-SP system. Panel (a) 
corresponds to the value of ^(Ra) in the parameter space 
{z, ACu} to distinguish the strong and weak coupling regions. 
Panel (b) and (d) show the dynamics of the excited state 
population and the two-photon correlation function {g^'^\t)) 
respectively for the three points plotted in panel (a). Panel (c) 
shows the SSQ luminescence spectra for those three particular 
cases. 



of the system not as a reservoir but as an ingredient co- 
herently coupled to one or more SSQs^^. 
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